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The vast majority of biochemical systems involve the exchange of information between different 
compartments, either in the form of transportation or via the intervention of membrane proteins 
which are able to transmit stimuli between bordering compartments. 

The correct quantitative handling of compartments is, therefore, extremely important when mod- 
elling real biochemical systems. The Bio-PEPA process algebra is equipped with the capability of 
explicitly defining quantitative information such as compartment volumes and membrane surface ar- 
eas. Furthermore, the recent development of the Bio-PEPA Eclipse Plug-in allows us to perform a 
correct stochastic simulation of multi-compartmental models. 

Here we present a Bio-PEPA compartmental model of the cAMP/PKA/MAPK pathway. We 
analyse the system using the Bio-PEPA Eclipse Plug-in and we show the correctness of our model by 
comparison with an existing ODE model. Furthermore, we perform computational experiments in 
order to investigate certain properties of the pathway. Specifically, we focus on the system response 
to the inhibition and strengthening of feedback loops and to the variation in the activity of key path- 
way reactions and we observe how these modifications affect the behaviour of the pathway. These 
experiments are useful to understand the control and regulatory mechanisms of the system. 



1 Introduction 

Compartments are widely present in biological systems and play a major role in a lot of biological 
processes [1]; they are enclosed by membranes, which isolate biological species from the external en- 
vironment. Biological species can be either located inside a compartment or reside on a membrane; 
in the latter case they can be attached to a side of the membrane (receptors, for instance) or span the 
entire membrane (i.e. transmembrane proteins). Species located in the same compartment can interact 
with each other through biochemical reactions; species can also move from one compartment to another 
via diffusion or through membrane channels, thus allowing the passing of information between adjacent 
compartments. Alternatively, the transmission of signals between compartments can involve receptors, 
which respond to the input of signalling molecules on one side of the membrane by triggering a cascade 
of events on the other side. 

Due to the fundamental importance of compartments, their correct handling is often essential when 
studying biological systems, and an increasing number of modelling languages support their explicit 
definition [29, 5, 33, 31, 6, 41]. These languages differ from each other in the kind of abstraction, the 
assumptions made (for instance, static or dynamic compartments) and the features they are designed to 
represent. Most of them are supported by analysis tools [40, 28, 39, 13]. 

In this work we consider Bio-PEPA with locations as our modelling language [11]. Bio-PEPA loca- 
tions abstract both biological compartments and membranes: they are defined by names, enriched with 
additional information allowing the modeller to express their kind (i.e. compartment or membrane), their 
size, and their position with respect to the other locations of the system. 
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Bio-PEPA offers a high level of abstraction, similar to that used to describe biochemical networks 
in the literature and in biological databases (see, for instance, models in the BioModels database [25]). 
Only static locations can be represented (i.e. compartments cannot merge, split, or undergo any structural 
change); nevertheless, their size can vary with respect to time. 

The main aim of this work is to demonstrate the power of Bio-PEPA with locations as a modelling 
language for multi-compartmental biochemical networks. Thus, we consider a moderately sized sig- 
nalling pathway describing the cAMP/PKA/MAPK activity [27] with the explicit definition of compart- 
ments and membranes and we translate it into Bio-PEPA. In this pathway the spatial aspects are partic- 
ularly relevant, since it involves interaction between molecular species residing in different biological 
compartments with different volume sizes. 

The analysis of the model is performed using the Bio-PEPA Eclipse Plug-in [10], a novel simulation 
tool for the construction and analysis of Bio-PEPA models. The Bio-PEPA Eclipse Plug-in offers meth- 
ods for numerical integration of Ordinary Differential Equations (ODEs) and stochastic simulation. We 
present here the latest version of the tool, which supports the definition of multi-compartmental models 
defined using the Bio-PEPA extension with locations and is able to handle multiple compartments with 
different sizes in the appropriate way. 

Using the Bio-PEPA Eclipse Plug-in, we validate the model against the expected results in the lit- 
erature and we perform in-silico experiments in order to investigate some properties of the pathway. 
Specifically, we study how the behaviour of the pathway is affected by varying the enzymatic activity of 
three key molecular species and by the removal of feedback loops with an important regulatory effect 
over particular species. These experiments are useful for a better understanding of the regulatory and 
control mechanisms of the system. 

The structure of the paper is as follows. Section 2 reports a brief introduction of the Bio-PEPA 
language, while Section 3 is devoted to the description of the tool. In Section 4 we present a model of the 
cAMP/PKA/MAPK signalling pathway, and in Section 5 we show the validation of the model and some 
analysis results. Section 6 concerns some related work, whereas the last section reports some conclusive 
remarks. The full Bio-PEPA model of the pathway is reported in Appendix A. 

2 Bio-PEPA with locations 

In this section we describe the main features of the extension of Bio-PEPA [12] with locations [11], the 
modelling language used in this work. 

A Bio-PEPA system representing a biochemical network consists of a context, defining information 
such as locations, kinetic rates, parameters and auxiliary information, a set of species components and a 
model component. 

Locations represent biological compartments and membranes. They are described by names, en- 
riched with additional information allowing the modeller to express their position with respect to the 
other locations of the system, their kind (i.e. compartment or membrane), and their size (i.e. volume or 
surface area). The structure of the biological system is modelled as a static hierarchy, represented as a 
tree whose nodes represent locations; each node has one child for each of its sub-locations. The location 
tree allows us to keep track of the relative positions of locations and must be associated with the location 
definition. 

The locations are defined as follows. 

Definition 1 Each location is described by "L : size unit, kind", where L is the (unique) location name, 
"size" expresses the size and can be either a positive real number or a more complex mathematical 
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expression depending on time t; the (optional) "unit" denotes the unit of measure associated with the 
location size, and "kind" e |M,C| expresses if it is a membrane or a compartment, respectively. The set 
of locations is denoted by X- 

Functional rates are used to describe kinetic laws associated with a reaction. They have the form 
"fja = expr", where a is the action type abstracting the reaction and expr is a mathematical expres- 
sion describing the law. The kinetic parameters used in the mathematical expressions are defined by 
parameters in Bio-PEPA. 

Species components describe the behaviour of the individual biological species, whereas the model 
component defines the interactions between the various species and the amount of them present in the 
system. The syntax of the Bio-PEPA components is defined as: 

5 ::= (ar,#f) op 5 | S +S | C \ S@location with op = J, | T I © I 9 I O P ::= S(x) \ P*fP 

where 5 is the species component and P is the model component. In the prefix term (a, k) op S , k is 
the stoichiometry coefficient of species S in reaction a, and the prefix combinator "op" represents the 
role of S in the reaction. Specifically, J, indicates a reactant, f a product, © an activator, 8 an inhibitor, 
and O a generic modifier. We can use "a op S" as an abbreviation for "(a, 1) op S". The operator "+" 

def 

expresses the choice between possible actions, and the constant C is defined by an equation C - S . The 
term S ©location indicates that the species represented by the process S is the location location. The 
parameter x e R + in S(x) represents the number of molecules, which can be abstracted to a level [9] for 
analysis as a CTMC. Finally, the process P^fQ denotes the cooperation between components: the set I 
determines those activities on which the operands are forced to synchronise. The notation >c means that 
the species components are obliged to synchronize on the common action types. 

Note that the same species in different locations is represented by distinct species components. There- 
fore, in order to avoid the possible duplication of actions in the model, we propose a notation to represent 
a species in multiple locations in a compact way. Given a species S which can be in L\, . . .,L n , one sin- 
gle component S is defined. This is simply a shorthand for a set of definitions S@L\,...,S@L n , each 
of which contains only those actions which can occur in the respective location. A reaction a\ occur- 
ring only in L; is defined as (<u\,k) op S @U, while a 2 occurring in all locations L\,. . .,L n is defined as 
(a 2 ,K) op 5. 

No constraint is imposed on the location of interacting species. This implies that it is possible to 
describe reactions involving species located in non-adjacent locations. The reason behind this choice 
is that, even though interactions are generally local (i.e. involve molecules in the same or adjacent lo- 
cations), imposing constraints on the relative position of species locations could prevent the modellers 
from defining some reactions (e.g. a reaction abstracting a sequence of reaction steps involving various 
locations). 

The formal definition of the Bio-PEPA system is the following. 

Definition 2 A Bio-PEPA system f is a 6-tuple {H,N ,"K,Tr, Comp, P), where: X is the set of locations, 
N is the set of (auxiliary) information for the species, 7C is the set of parameters, fj{ is the set of 
functional rates, Comp is the set of species components and P is the model component. 

A main feature of Bio-PEPA is that it can be seen as a formal, intermediate representation of bio- 
chemical systems, on which different kinds of analysis can be carried out, through the defined mappings 
into continuous-deterministic and discrete- stochastic modelling languages and CTMC with levels. 

Bio-PEPA is given an operational semantics, which is used for the derivation of the stochastic la- 
belled transition system and the associated CTMC with levels. We do not report here the details on 
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the operational semantics since they are not needed for the kinds of analysis on which we focus in this 
work (continuous-deterministic and discrete-stochastic). For details, see [11, 12]. 

Finally, it is important to remember that Bio-PEPA (and, in general, process algebras) offers a com- 
positional approach for the definition of models: it is possible to build the full model starting from the 
definition of its subcomponents. This approach is particularly useful for biological systems, as they are 
generally large and complex. For these systems, a compositional definition of the model is desirable 
since it makes the construction of the full model easier and it helps in understanding the relationships be- 
tween its various subcomponents and the structure of the network. The compositional approach has been 
used in the definition of the Bio-PEPA model for the cAMP/PKA/MAPK pathway. The compositional 
rules for Bio-PEPA models are standard and not reported here. 

3 The Bio-PEPA Eclipse Plug-in 

The Bio-PEPA Eclipse Plug-in [2] is one of two tools developed at Edinburgh University for the Bio- 
PEPA language. Built on top of the Eclipse platform [16], the plug-in offers an operating system agnostic 
editor and suite of time-series solvers for the construction and analysis of Bio-PEPA models. It provides 
the modeller with a rich modelling environment for Bio-PEPA, from the initial creation of the model 
through to the analysis and inspection of the results. 

The features currently supported in the plug-in can be roughly grouped into those assisting the mod- 
eller in the writing of their models and features for the analysis of the model. Examples of the former 
include an editor (with Bio-PEPA-aware syntax highlighting) and static analysis of the model for detect- 
ing various syntactic errors. For the latter, the plug-in is capable of generating and displaying time-series 
results, from ODE numerical integration or stochastic simulation algorithms, with the option to export 
these results in both graphical and textual formats. Many of these features can be seen in Figure 1. 

As with all previous releases of the plug-in, the language is augmented with keywords for particular 
definitions. These keywords increase the readability of the model, clearly signposting those definitions 
pertaining to locations, species and kinetic laws. These indicators are then emphasised further by the 
syntax highlighting built into the developing environment. Together they assist the modeller in recog- 
nising the purpose of various definitions found within any given model. As of version 0. 1 .0, the Eclipse 
Plug-in will contain initial support for multi-compartmental models as defined by the Bio-PEPA exten- 
sion with locations. The following description details both this support along with several refinements 
to the implemented language. For the sake of brevity the description does not cover the entire language, 
instead concentrating on those required to model the system of interest. 

Units of measures are currently unsupported, all values must be defined with respect to the same 
generic unit, with compartments referring to volumes (unit 3 ) and membranes as surface area (unit 2 ). Fi- 
nally, the initial values for species are assumed to be in terms of molecular counts. With the introduction 
of locations with different sizes, concentrations become an inappropriate form of measurement and must 
be translated into molecular counts. 

3.1 Supported Language 

Bio-PEPA does not impose a strict ordering on the definitions of context and species components, but 
these must all precede the model component. This final statement provides the structure of the system as 
well as the initial conditions. Since the modeller is relatively unconstrained with respect to the ordering 
we shall present the syntax for different definitions in the order that they are parsed within the plug-in. 
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Figure 1 : Screen shot of the Bio-PEPA Eclipse Plug-in. The editor can be seen in the top left corner, with 
the outline view in the top right. The problems tab (bottom right) displays warnings and errors associated 
with the open models and the ability to plot results can be seen in the bottom left. 



The first is the location definition, an example of which can be seen below: 

location cyto in cytojnem : size = 1, kind - compartment; 
location cytojnem : size = 0.2, kind = membrane; 

The required keyword location identifies this as a location statement and the optional keyword in 
is used for defining the position of the given location in the location tree. In the example above we 
have defined a compartment cyto to exist within a second location known as cytojnem. Allowable 
parameters for these locations include the size and kind, both of which can be seen. For reasons of 
legibility acceptable kinds are compartment and membrane rather than M and C. Based on the kinds 
in the example and the spatial hierarchy we can clearly see that cytojnem is the membrane to the cyto 
compartment. 

The next required definition in a Bio-PEPA model describes the kinetic rates of the reactions, with 
an example below: 

kineticLawOf vlf : kf jactivateJ3s*isoJBARjG@ cytojnem; 

where kf .activate JGs is the constant rate of the reaction vlf and isoJ}ARjG@cytojnem is the reactant 
of the reaction (the species amount is given in number of molecules). Like location definitions, this def- 
inition is easily identifiable by the required keyword kineticLawOf, followed by its label for reference 
within the species component definitions. The defined rate can either be a custom rate, as seen in the 
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example, or make use of predefined laws such as mass-action or Michaelis-Menten. For example, if 
we know that the only reactant in the reaction associated with vlf is iso_BARJ3@cytojnem, then this 
custom rate describes mass-action. This could then be re-written as: 

kineticLawOf vlf : fMA(kf .activate. Gs); 

where fMA(r) stands for mass-action kinetic law with constant parameter r. The inclusion of predefined 
laws offers two clear benefits to the modeller. Mistakes in inputting rates can be reduced by removing 
the need to input the laws manually. Secondly, the use of predefined kinetic laws allows for additional 
static analysis of the model. If a reaction is defined as using Michaelis-Menten (fMM(v m,^m), where 
vm, kM are the two constants in Michaelis-Menten kinetic law 1 ), then the reaction must have a single 
species acting as a reactant or substrate, another species performing the role of the catalyst or enzyme 
and a third as the product of the reaction. If the list of species involved in this reaction differs from those 
expected for the kinetic law stated the tool will flag this as an error for the modeller. 

For the species component definitions, acceptable symbols for the different possible behaviours are 
required. Below is the list of Bio-PEPA operators along with their ASCII equivalent for creating Bio- 
PEPA models. 



Behaviour 


Bio-PEPA symbol 


ASCII representation 


reactant 




« 


product 


t 


» 


activator 


© 


(+) 


inhibitor 


e 


(-) 


modifier 


o 


(•) 


uni-directional transportation 


— > 


- > 


bi-directional transportation 


<-> 


< - > 


cooperation 


X 

I 





These operators are used within the species component definitions, an example of which can be seen 
below. 

PDE4 = v27 » PDE4@cyto + v20 « PDE4@cyto + v5 (+)PDE4@cyto; 

The species PDE4 is involved in three interactions with different roles: it is a product in the reaction 
v27, a reactant in v20 and an enzyme in v5. In the example above all three reactions are only possible 
when the species PDE4 is in the location cyto. From this example we can see shorthand notations for 
when the stoichiometric coefficient is one. As with all defined shorthands, their purpose is to keep the 
definitions free of superfluous details. This is not to suggest that stoichiometric information in general is 
unnecessary, but that through the use of this shorthand we can emphasise the less common instances (a 
stoichiometry of one being the majority of observed values) where this is not the case. 

The final definition in a Bio-PEPA system is that of the model component (see Section 2). The 
model component describes the initial amount of each species in the system and the possible interactions 
involving these species. The plug-in however offers true compositionality. It does this by first augmenting 
the model component syntax to accept names and also allowing definitions written in the syntax of the 

1 Given an enzymatic reaction S +£ — > P + E, where 5 is the substrate, P is the product and E is the enzyme, the Michaelis- 
Menten kinetic law is vm • E ■ S / (Km + S ), where vm is the maximum rate constant and Km is the Michaelis constant. 
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model component to be assigned names. These additional labelled compositional definitions can then be 
referred to in either the model component or other labelled compositional definitions. This allows the 
modeller to describe a complex system in terms of smaller pathways, or indeed any logical grouping that 
they wish to enforce. 

MAPK_Pathway ::= B_Raf@cyto[\20] <*> B_Rafjactive@cyto[0] <*>... 

<*> MAPK@cyto[2l7] <*> MAPK_active@cyto[0] 
<*>•••<*> PTP_PP@cyto[60]; 

Here we have a labelled composition definition for the MAPK pathway. It should be noted that la- 
belled compositional definitions use a different assignment operator to other definitions, using :: = instead 
of =. As stated earlier, and as can be seen in the example above, labelled compositions also use the defi- 
nition terminator. This leaves the model component as the final definition to present, which an example 
of can be seen below: 

G_ Pathway < v09,vl7 > c AMP .Pathway < v08,vl5 > MAPK .Pathway 

The benefits gained through compositional definitions are highlighted in this model component def- 
inition. Not only does the model component make it clear that the system is compromised of three 
pathways, but can also serve to highlight important reactions. We can see that interaction between the 
cAMP and MAPK pathways is limited to reactions v08 and vl5. If the MAPK pathway is reliant on the 
cAMP pathway, rather than any interaction being redundancy in the system, then the rates for these reac- 
tions are likely candidates for perturbation. This of course is not guaranteed, merely that the cooperation 
sets can be used to highlight limited interactions between groups of species. The compositional approach 
also makes it very easy to enable or disable entire pathways. 

Once the model has been built, the plug-in parses the model and performs static analysis to detect a 
variety of modelling errors, an example of which was given earlier in regards to kinetic rate laws. Once 
the model is free of these errors the modeller has access to ODE solvers in the form of a Runge-Kutta 
implicit-explicit solver and a Dormand-Prince adaptive step-size solver [38], these being made available 
through the odeToJava library. The plug-in also supports analysis through stochastic simulators, with 
implementations of Gillespie's Stochastic Simulation Algorithm (SSA) [19], Gibson-Bruck [17] and 
the Tau-Leap algorithms [18] available. The ability to use different types of solvers has allowed us to 
discover errors in published modelling studies in computational biology [4]. 

4 The cAMP/PKA/MAPK pathway 

The cAMP/PKA/MAPK network is a ubiquitous signalling pathway. In this work we consider the system 
as modelled in [27] and in particular the ODE model described by the authors and available in SBML 
format from [3]. 

The cAMP-dependent protein kinase (PKA) is responsible for performing several functions in the 
cell, such as regulation of glycogen, sugar, and lipid metabolism. Among the downstream targets of 
PKA there are mitogen-activated protein kinases (MAPKs), which are essential in regulating various 
cellular activities, including gene expression, mitosis, differentiation, and cell survival/apoptosis. 

The major role of this pathway in various cellular processes, together with the availability of quan- 
titative data such as the sizes of the involved compartments, make it an interesting case study for our 
purposes. 
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As is common in signalling pathways, spatial aspects are particularly relevant in this system, since 
it involves interaction between molecular species residing in different biological compartments. We 
consider molecules lying inside the cell (i.e. in the cytoplasm), and outside it (i.e. in the extra-cellular 
environment), specifically the receptor agonist which triggers the signalling cascade. These two com- 
partments are separated by a membrane (i.e. the cytoplasmic membrane) where receptors and specific 
membrane-proteins are located. 

The compartments have different sizes and this exercises a direct impact on the quantitative behaviour 
of the system. As in the original ODE model [27, 3], we consider the cell to be a sphere with a surface- 
to-volume ratio equal to 0.2. The kinetic parameters we use are taken from the original paper, which 
have been derived from the biochemical literature. 

Before presenting the Bio-PEPA model in the following Section 4.1, we briefly describe the biologi- 
cal system. Figure 2 shows a schematic representation of the pathway, and the full model is reported in 
Appendix A. 
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Figure 2: Schematic representation of the pathway. The coloured boxes represent a decomposition of the 
pathway into three subnetworks: the G-protein mediated activation of AC, the activation of PKA, and 
the MAPK signalling cascade. Arrows represent reactions, and lines terminating with circles represent 
enzymes participating in reactions. 

The cAMP/PKA/MAPK network we consider is in neurons, and its biological action is initiated 
by isoproterenol (ISO) binding to specific cell-surface receptors called /J-adrenergic receptors (BARs), 
which are a class of G-protein coupled receptors (GPCRs). 

GPCRs reside in the cytoplasmic membrane and respond to a wide range of extracellular stimuli, 
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in our case to the synthetic drug isoproterenol. When GPCRs bind an extracellular ligand, they interact 
with G-proteins and propagate the signal across the membrane into the cytoplasm. 

G-proteins are heterotrimeric complexes bound to the inside surface of the cell membrane, and con- 
sist of an active subunit G a and two tightly associated Gp and G y subunits. G a subunits alternate between 
an inactive guanosine diphosphate (GDP) and an active guanosine triphosphate (GTP) bound state. The 
transition from the inactive to the active state of G a and its simultaneous dissociation from the G y p 
subunit and the GPCR are stimulated by the binding of GPCRs with the ligand. 

G s subunits are a subclass of G a subunits, and they stimulate the production of cyclic-AMP (cAMP) 
from ATP, by activating (by binding to it) the membrane-associated enzyme adenylate cyclase (AC) 
which, in turn, produces cAMP. 

cAMP is a key signalling molecule which acts as a second messenger, activating the cAMP-dependent 
protein kinase A (PKA). PKA is a tetrameric holoenzyme, consisting of two regulatory and two catalytic 
subunits (/?2C2). Under low levels of cAMP, PKA is inactive; when the concentration of cAMP rises, 
cAMP molecules bind to each of the two binding sites on the two regulatory subunits, which leads to 
the dissociation and activation of the catalytic subunits (i.e. four cAMP molecules are needed to activate 
one PKA molecule). Upon activation, the catalytic subunits can phosphorylate a number of downstream 
targets. 

PKA is an enzyme which participates in several metabolic and signalling pathways, leading to a wide 
range of biological responses. Here we focus on its ability to initiate the pathway leading the activation 
of the Mitogen- Activated Protein Kinase (MAPK). 

The downstream pathway in our network represents an instance of the common MAPK signalling 
chain. It starts with the phosphorylation of the MAPKKK B-Raf, which phosphorylates (i.e. activates) 
the MAPKK MEK, which in turn phosphorylates MAPK. 

The amounts of PKA and MAPK in cells are important elements to monitor when analysing the 
behaviour of the system, because of their essential role in several signalling pathways and other cellular 
processes. 

An interesting aspect of this pathway is the presence of numerous feedback mechanisms, which are 
carried out mainly by kinases and phosphatases, regulating the amount of the involved proteins. 

Two kinds of phosphatases are involved: first, the serine/threonine protein phosphatase 2 (PP-2A), 
which is a ubiquitous enzyme able to dephosphorylate a wide range of proteins among which there are 
the components of the MAPK signalling cascade (B-Raf, MEK and MAPK); second, protein tyrosine 
phosphatases (PTPs), which are also involved in the dephosphorylation of MAPK. 

In addition to its positive role carried out by activating B-Raf, PKA is further involved in the reg- 
ulation of the pathway, by exercising both a positive and a negative regulatory role: respectively, it 
phosphorylates PTPs (in the phosphorylated form, PTPs are less effective in dephosphorylating MAPK) 
and PDE4 (in the phosphorylated form, PDE4 enhances the transformation of cAMP into AMP, making 
it unavailable for activating PKA). 

The relative strength of the different PKA enzymatic activities, and particularly the opposed positive 
and negative feedbacks, is therefore an interesting issue to be considered when studying PKA regulatory 
role. 

4.1 A Bio-PEPA Model of the Network 

In the schematic representation of the pathway in Figure 2 we have subdivided the system into three 
different subnetworks: 

• upstream pathway: G-protein mediated activation of AC 
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• intermediate pathway: activation of cAMP-dependent PKA 

• downstream pathway: MAPK signalling cascade 

The Bio-PEPA language allows the definition of compositional models and, therefore, the above sub- 
division can be straightforwardly employed in the model itself. This helps in improving the readability 
of the model. The full Bio-PEPA model (using the syntax and the shorthand notation supported by the 
Bio-PEPA Eclipse plug-in) is reported in Appendix A. 

Three locations are defined in the Bio-PEPA model: cyto, cytojnem, and extra. Each of them 
represents one of the involved biological compartments, and is assigned a size expressing its volume (3D 
compartments) or its surface area (2D membranes). 

Kinetic parameters and rate laws are specified. Some of the reaction rates follow the mass-action 
law, while others are based on the Michaelis-Menten kinetics. The quantitative data we adopt are those 
defined in the original ODE model [3]. The reversible reactions are decomposed into two irreversible 
(mass-action) reactions describing the forward and the backward directions, respectively. The full set of 
parameters used is reported in the supplementary material. 

Each molecular species in the pathway is abstracted by a species component describing the reactions 
in which it is involved and its role in each of them (i.e. whether it is a reactant, a product, or an enzyme). 
The Bio-PEPA model component describes how the species components interact with each other, and 
defines the initial amounts for each species. 

There are two noteworthy differences with respect to the original model. These are related to the 
conversions needed to switch to molecular counts as units of measure for the species amount. In the 
ODE model the species amounts are specified in terms of concentrations (numbers of molecules over 
volume/area), while in our Bio-PEPA model they must be specified in terms of numbers of molecules. 
Consequently, the rate laws must be rescaled accordingly to take into account the sizes of the compart- 
ments in which the species involved in the reactions are located. 

In this definition of the model, compartment volumes and membrane surface areas are assumed to 
be constant. This constraint is not present in Bio-PEPA, which allows location sizes to be defined as 
arbitrary functions (e.g., they can change over time or depend on the amount of given species). Here 
we do not consider generic location sizes in order to be able to compare our results with the original 
model. With this assumption, the system could be modelled using languages which are not equipped 
with a notion of location (provided reaction rates are carefully computed taking into account volumes). 
However, using a language, such as Bio-PEPA, which is equipped with a notion of location, has also the 
obvious advantage that the explicit definition of information such as locations of species, volumes and 
areas brings a greater clarity of the model. 

5 Model Analysis 
5.1 Validation 

In order to verify that the Bio-PEPA model is a faithful representation of the original ODE model [3], we 
compared the results we have obtained from the Bio-PEPA model with the results obtained through the 
BioModels online simulator on the original model. 

Figure 3(a) reports the comparison of the results for the ratio of active MAPK (i.e. phosphorylated 
over total MAPK). We show both the time series obtained using one of the stochastic simulators and the 
one obtained using one of the ODE solvers available within the Bio-PEPA Eclipse Plug-in. We observe 
a good agreement between the results from the Bio-PEPA model and the ones from the original ODE 
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model. In particular, the two ODE-based analyses produce coincident results. The SSA plot, obtained as 
the average over 500 individual stochastic simulation runs, is in good agreement with the deterministic 
analysis and also shows the variability of the system due to the presence of some stochastic noise. 



MAPK-P / (MAPK + MAPK-P) 




Intermediate products 










— AC active 




— PKA active 




— MAPK active 











100 200 300 400 500 600 700 
Time 



(a) Ratio of active MAPK 



(b) Active AC, PKA and MAPK 



Figure 3: Comparison between the simulation results obtained by the Bio-PEPA Eclipse Plug-in and the 
results of the original ODE model. SSA results are the average over 500 simulation runs. 

A similar agreement was also found for all the other molecular species involved in the pathway. In 
Figure 3(b) we report the results for three of the key molecular species involved in the pathway (AC, 
PKA, and MAPK), each of which is the final product of one of the three subnetworks. For each of these 
proteins, we compare the amount of its active form in our SSA results (solid lines) with the results from 
the original model (dashed lines). 

Note that, in the SBML model, the amounts of species are given in terms of numbers of molecules 
over compartment sizes and, hence, the amounts must be rescaled by multiplying by the compartment 
sizes in order to be comparable with our results, which are expressed in terms of molecular counts. AC is 
located on the cell membrane and, therefore, it must be multiplied by the area of the membrane (0.2 jura 2 ); 
in the case of PKA and MAPK the amounts result the same, since the volume of the cytoplasm is 1 yum 3 . 



5.2 Experimentation 

In this section we report some of the computational experiments we have performed on the model. First, 
we have studied the effect, on the final product of the pathway (MAPK-P, i.e. phosphorylated MAPK), 
of changes in the network structure: we have removed some of the feedback loops, in particular the ones 
involving PKA, which is one of the proteins mostly responsible for both positive and negative regulation 
of various pathway components. 

Figure 4 shows the effect of inhibiting some important steps on the behaviour of the system. In Fig- 
ure 4(a) we compare the behaviour of the original system (blue line) with the behaviour of three modified 
systems: one with weaker enzymatic activity of PKA (10-fold increase of K m _v0&) in the phosphoryla- 
tion of B-Raf (red line), one with slower activation of AC (10-fold decrease of K f _AC .activation in 
the binding reaction between AC and G as : GTP) (green line), and one with slower cAMP/R2C2 bind- 
ing (10-fold decrease of Kf_vl3, Kf_vl&, Kf_v23, Kf_v32) (black line). Each of these modified systems 
corresponds to slowing down one reaction rate in each of the subnetworks; Figure 4(a) shows that a re- 
duction in the rate of cAMP binding has the strongest effect: an expected result, since the cAMP binding 
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reaction must occur four times for a single molecule of PKA to become active. 
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Figure 4: Experimenting on reaction rates and feedback loops (each line is the average behaviour over 
500 independent stochastic simulation runs). 

In Figure 4(b), instead, we experiment with the model by changing the parameters relative to the 
enzymatic activity of PKA and studying their effect on the system behaviour. As described in Section 4, 
PKA takes part in a positive feedback loop by phosphorylating PTPs, and in a negative feedback loop by 
phosphorylating PDE4. 

The results shown in Figure 4(b) confirm the presence of these feedback loops: a weaker enzymatic 
activity of PKA (100-fold increase of K m _vl5) in the phosphorylation of PTP (red line) leads to a de- 
crease of MAPK-P (due to the lower presence of phosphorylated PTP, which is less effective than its 
unphosphorylated form in dephosphorylating MAPK), while a stronger enzymatic activity of PKA (100- 
fold decrease of ^T m _vl5) (green line) leads to an increase of MAPK-P. The opposite occurs for the 
phosphorylation of PDE4: a weaker enzymatic activity of PKA (100-fold increase of K m jv20) in the 
phosphorylation of PDE4 (black line) leads to an increase of MAPK-P (due to the lower presence of 
phosphorylated PDE4, which is more effective than its unphosphorylated form in transforming cAMP 
into AMP), while a stronger enzymatic activity of PKA (100-fold decrease of K m _v20) (cyan line) leads 
to a decrease (though smaller) of MAPK-P. 



6 Related work 

Several languages have been defined to model biological systems with compartments and membranes [29, 
5, 33, 31, 6, 41, 24]. In the following we report a brief overview of the most well-known languages and 
the associated analysis tools. 

Membrane systems or P systems [29, 8] are computational models that are based on a notion of 
membrane structure. In particular, multisets of objects are enclosed in a nested hierarchy of membranes 
and their behaviour is described by local rewriting rules. 

The analysis of P systems is supported by the P System Modelling Framework [28]. This framework 
provides stochastic simulation based on the multi-compartmental Gillespie algorithm [30] and offers the 
mapping to PRISM [22, 32] for model checking. Furthermore, translations to SBML [35] and CellDe- 
signer [7] have been implemented in order to export P system models into other formats for further kinds 
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of analysis and for graphical representation. 

The variant of membrane systems with peripheral and integral membrane proteins defined by Cava- 
liere and Sedwards [6] is at the basis of Cyto-Sim [36, 13]. The analysis is limited to stochastic simulation 
by Gillespie's algorithm [19], but it is possible to export the corresponding ODE model into MATLAB 
format [26] for further analyses. In addition to mass-action kinetics, the simulator supports chemical 
reactions with arbitrary kinetic laws based on functions of the reactants. 

One of the first process calculi with an explicit notion of compartments is BioAmbients calculus [33]. 
A BioAmbient system is seen as a hierarchy of nested ambients, generally abstracting compartments, 
containing communicating processes whose actions specify the evolution of the system. Various kinds 
of action involving compartments can be easily represented, such as transport of small molecules across 
compartments. BioAmbients models are supported by two tools: Bio-SPI [37], a stochastic simulator 
based on Gillespie's algorithm and originally implemented for models in stochastic /r-calculus, and the 
BioAmbient machine (BAM) [43, 40], a more efficient and user-friendly simulator written in Java. BAM 
is based on the stochastic version of BioAmbients, defined by Vigliotti and Harrison [42] and can be used 
to simulate bio-regulatory pathways as well as membrane interaction. BAM automatically produces both 
a simulation graph and a debugging trace of the program. Furthermore, it offers a graphical user interface 
for entering program code and visualising simulation results. 

A more active role for membranes can be described using Brane calculi [5]. In Brane calculi, a 
system is represented as a set of nested membranes, and a membrane is represented as a set of actions. 
Membranes can move, merge, split, enter into and exit from other membranes. 

The bioK-calculus [24] is a language for describing proteins and cells based on Brane calculi and 
/c-calculus [14]. The atomic elements of the calculus are proteins and two constructors for representing 
solutions and cells. Proteins are characterised by sites, that can be in different states: bound to another 
site of a protein, visible (i.e not connected to other sites) and hidden (i.e not connected to other sites 
but not available for other interactions). Cells abstract compartments and consist of a membrane and a 
cytoplasm. Reactions can be interactions between two proteins or between two membranes. 

In Beta-binders [31] and the associated BlenX language [15], systems are modelled as a composition 
of boxes representing biological entities. In the original version of the language compartments are not 
defined explicitly, but a virtual form of nesting is rendered by appropriate typing for sites. Explicit static 
compartments and transport of biological across them have been added to Beta-binders in [21]. The 
analysis of BlenX models is possible in the Beta Workbench (BetaWB) [39]. The BetaWB is composed 
of a stochastic simulator based on an efficient variant of Gillespie's algorithm, the BetaWB designer, a 
graphical editor for developing models and the BetaWB plotter, a tool to analyse the results of a stochastic 
simulation run. Furthermore, BlenX models can be exported in SBML. 

Finally, there are some variants of the stochastic /r-calculus enriched with the notion of locations 
(see [41, 23], for instance). S/r@ calculus [41] extends the stochastic /r-calculus syntax with the ex- 
plicit addition of compartments. This language handles varying volumes and dynamic compartments by 
defining the compartment volume as the sum of the volumes occupied by all the molecules it contains. 

Another extension of the /r-calculus with compartments is SpacePi [23]. Pi processes are embedded 
into a vector space and move individually. Only processes that are sufficiently close can communicate. 
The operational semantics of SpacePi defines the interplay between movement, communication, and 
time-triggered events. 

All the languages listed above differ from Bio-PEPA in various aspects: they are based on differ- 
ent levels of abstraction or focus on dynamical compartments or handle volumes in a distinct way. In 
Bio-PEPA, locations have a fixed structure and compartments are essentially containers for biological 
species: interactions between molecules in different compartments is allowed, but the main evolution is 
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given by interactions of molecular species within compartments. This concept of location is that usu- 
ally considered in biochemical networks in databases and in the literature. For these systems Bio-PEPA 
offers a direct and formal representation, allowing an intuitive representation of both intra-compartment 
and inter-compartment reactions. 

Concerning the Bio-PEPA Eclipse Plug-in, a prominent feature of the tool is its support of both 
stochastic simulators and ODE numerical solvers. The ability to utilise different analysis techniques is 
very important; this can aid in understanding different behavioural aspects of the system and discovering 
possible errors due to the use of a specific solver/simulator [4]. Furthermore, the user can select the most 
appropriate analysis method based on the model under study. Most of the other tools support stochastic 
simulation by Gillespie's algorithm only and in a few cases mappings to other languages or tools, such 
as SBML or PRISM. 

7 Conclusions 

Interactions of molecular species located in compartments with different sizes occur typically in biolog- 
ical systems. In signalling pathways, for instance, a biochemical signal is transferred from the extra- 
cellular environment into the cell through the intervention of membrane receptors and other membrane- 
bound proteins. 

In this work we have presented a Bio-PEPA compartmental model of the cAMP/PKA/MAPK path- 
way. Our main aim was to demonstrate, on a real case study, the effectiveness and the correctness of the 
representation of multi-compartment models in Bio-PEPA. 

The Bio-PEPA language allows us to explicitly define locations representing compartments and 
membranes, to express the location of each species and to indicate the relative position of a location 
with respect to the others. Furthermore, Bio-PEPA correctly deals with the location sizes (volumes and 
surface areas). It is generally useful to collect these kinds of information, in order to have a better view 
of the system under consideration. Moreover, the information about the size of locations is necessary to 
compare the results obtained in Bio-PEPA with the ones in the literature and experiments, when given in 
concentrations instead of number of molecules. Another benefit is although the system in question can 
be written as a set of ODEs, high-level languages like Bio-PEPA can offer a cleaner view of the model 
through a more restricted (but focussed) number of operators. Through these operators the modeller must 
explicitly state any behaviours he or she wishes the components of the system to exhibit. A restricted set 
of operators also allows the tool to perform more extensive checking over any legal Bio-PEPA model, 
something not possible if manually coding a set of ODEs in a general tool such as MATLAB. 

In order to be able to validate our results, we developed our Bio-PEPA model as a faithful represen- 
tation of an existing ODE model. 

The analysis of the pathway was performed using the stochastic simulator available within the Bio- 
PEPA Eclipse Plug-in, which is able to handle multi-compartmental models. The comparison of our 
results with the ones obtained from the original ODE model demonstrate the correctness of our ap- 
proach. We also perform a number of simulation experiments in order to investigate specific properties 
of the system such as its response to the inhibition of feedback mechanisms and to the variation of key 
molecular species. 

In this model we have found it convenient to describe the behaviour of the pathway in terms of three 
functional elements, a feature supported by the tool. Currently such modularisation is supported only 
at the level of model definition as a convenience for the modeller. However it is an interesting area for 
future work to explore the extent to which such structuring might be exploited during model analysis. 
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While the locations extension for Bio-PEPA does not restrict species interaction by location, future 
improvements within the tool will flag such reactions. Once implemented this will allow the modeller 
to clearly see which reactions involve species not existing in adjacent locations. Other future work 
in the Bio-PEPA Eclipse Plug-in includes improved support for locations as well as the implementa- 
tion for automatic import/export of existing specification languages supporting multi-compartment mod- 
els (e.g. [35, 34, 20]). 

Finally, we would like to extend the analysis of this model (and of multi-compartmental models 
in general) by applying some of the CTMC -based analysis methods (e.g. model-checking) which are 
supported by the Bio-PEPA language. 
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A Bio-PEPA Model 

In this appendix we report the full Bio-PEPA model, specified in the tool syntax, corresponding to the 
cAMP/PKA/B-Raf/MAPKl,2 pathway. The abbreviations described in the paper are used. The measure 
unit for the volume size is pirn 3 and the unit for the area is pm 2 . The units of parameters are defined 
accordingly. 

In the definition of parameters, the constants omegaxyto and omega-extra are the scaling fac- 
tors used for converting parameters involving concentrations (expressed in pM) in the correct unit in 
terms of number of molecules, omegaxyto is for parameters of reactions in the cytoplasm whereas 
omega_extra is for parameters of reactions in the extracellular space, omega.cyto is defined as 10~ 6 -Na • 
size{cyto), where 10 -6 is to convert p M into M, Na = 6.022 • 10 23 is the Avogadro number (i.e. the num- 
ber of molecules in a mole of substance) and size{cyto) is the volume size of the cytoplasm. Similarly, 
omega_extra - 10~ 6 • Na • size{extra), where sizeiextra) is the volume size for the extracellular space. 

location cyto in cyto_mem : size = 1, kind = compartment; 
location extra : size = 0.111, kind = compartment; 

location cyto_mem in extra : size = 0.2, kind = membrane; 



omega_cyto = 602; 

kcat_pde4_p_pde4_p = 20; 

kcat_PDE4_PDE4 = 8; 

kcat_PKA_activates_Raf = 10; 

kcat_highKM_PDE = 8; 

kcat^C_basaLAC_basal = 0.2; 

kcat_PKA_P_PDE = 10; 

kcat_PTP_PKA = 0.1; 

kcat_PPase_MAPK = 0.636; 

kcat_pp_ptp_pp_ptp = 5; 

kcat_PPase_mek = 5; 

Kf_activate_Gs = 0.025; 

Km_pde4_p = 1.3; 

Kf_v4 = 1; 

Km_PDE4 = 1.3; 

KrbgbindsGRK = 0.5; 

Km_v08 = 0.5; 

Kf_GTPase = 0.067; 

Kf_trimer = 6; 

Kf_G_binds_iso_BAR = 10; 

Kf_vl3 = 8.35; 



omega.extra = 66.82; 

kcat JPPase Raf = 5; 

kcat_MEK_activatesJVIAPK = 0.15; 

kcat_AC_active_AC_active = 8.5; 

kcat_PKA_P_PTP = 0.2; 

kcat_grk_GRK = 0.104; 

kcatJRaLactivatesJVIEK = 0.105; 

kcat_PTP = 1.06; 

kcat_pp2a_4_pp2a_4 = 5; 

kcat_GRK_bg_GRK_bg = 1.34; 

Reactivate _Gs = 0; 

Km_v3 = 15.7; 

Kr_v4 = 0.2; 

Kf_bg_binds_GRK = 1; 

Km_v07 = 0.046; 

RmJVC_active = 32; 

Kr_GTPase = 0; 

Kr.trimer = 0; 

Kr_G_binds_iso_BAR = 0.1; 

Kr_vl3 = 0.017; 
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Km_vl4 

Kf_AC_activation 

Km_AC_basal 

Kr_vl8 

Km_v20 

Km_v22 

Kr_v23 

Kf_v25 

Km_v26 

Km_v28 

Kf_G_binds_BAR 

Km_v31 

Kr_v32 



= 


15; 


Km_vl5 


= 


0.1; 


= 


500; 


Kr _AC .activation 


= 


1; 


= 


1030; 


Kf_vl8 


= 


0.006; 


= 


0.00028; 


Km_grk 


= 


15; 


= 


0.5; 


Km_v21 


= 


0.159; 


= 


9; 


Kf_v23 


= 


0.006; 


= 


2.8e-4; 


Km_v24 


= 


0.46; 




1; 


Kr_v25 




0.062; 




0.77; 


Km_pp2a_4 




8; 




6; 


Km_GRK_bg 




4; 




0.3; 


Kr _G -binds _BAR 




0.1; 




15.7; 


Kf_v32 




8.350; 




0.017; 









kineticLawOf vOlf 
kineticLawOf vOlb 
kineticLawOf v02 
kineticLawOf v03 
kineticLawOf v04f 
kineticLawOf v04b 
kineticLawOf v05 
kineticLawOf v06f 
kineticLawOf v06b 
kineticLawOf v07 
kineticLawOf v08 
kineticLawOf v09 
kineticLawOf vlOf 
kineticLawOf vlOb 
kineticLawOf vl If 
kineticLawOf vl lb 
kineticLawOf vl2f 
kineticLawOf vl2b 
kineticLawOf vl3f 
kineticLawOf vl3b 
kineticLawOf vl4 
kineticLawOf vl5 
kineticLawOf vl6f 
kineticLawOf vl6b 
kineticLawOf vl7 
kineticLawOf vl8f 
kineticLawOf vl8b 
kineticLawOf vl9 
kineticLawOf v20 
kineticLawOf v21 
kineticLawOf v22 
kineticLawOf v23f 
kineticLawOf v23b 
kineticLawOf v24 
kineticLawOf v25f 
kineticLawOf v25b 
kineticLawOf v26 



fMA(Kf .activate _Gs) ; 

fMA(Kr_activate_Gs /(omega_cyto * omega_cyto)); 
fMM(kcat_pde4_p_pde4_p, Km_pde4_p * omega.cyto); 
fMM(kcat_PPase_Raf, Km_v3 * omega_cyto); 
fMA(Kf_v4 /omega_extra); 
fMA(Kr_v4); 

fMM(kcat_PDE4_PDE4, Km_PDE4 * omega_cyto); 
fMA(Kf_bg_binds_GRK / omega.cyto); 
fMA(Kr_bg_binds_GRK); 

fMM(kcatJVlEK_activates JVlAPK, Km_v07 * omegaxyto); 

fMM(kcat_PKA_activates_Raf, Km_v08 * omega_cyto); 

fMM(kcat_AC .active _AC .active, Km_AC_active * omega.cyto); 

fMA(Kf_GTPase); 

fMA(Kr_GTPase); 

fMA(Kf_trimer / omega_cyto); 

fMA(Kr_trimer); 

fMA(Kf_G_binds_iso_BAR/omega_cyto); 
fMA(Kr_G_binds_iso^AR); 
fMA(Kf_vl3 / omega_cyto); 
fMA(Kr_vl3); 

fMM(kcatJiighKM_PDE, Km_vl4 * omega_cyto) ; 
fMM(kcat_PKA_P_PTP, Km_vl5 * omega_cyto); 
fM A(Kf _AC _activation/omega_cy to) ; 
fM A(Kr _AC .activation) ; 

fMM(kcat_AC_basal_AC_basal, Km_AC_basal*omega_cyto); 

fMA(KLvl8 /omega_cyto); 

fMA(Kr_vl8); 

fMM(kcat_grk_GRK * cyto_mem / omega_cyto, Km_grk * cyto_mem); 
fMM(kcat_PKA_P_PDE, Km_v20*omega_cyto); 
fMM(kcat_Raf_activatesJVIEK, Km_v21 * omega.cyto); 
fMM(kcat_PTP_PKA, Km_v22 * omega.cyto); 
fMA(Kf_v23 /omega_cyto) ; 
fMA(Kr_v23); 

fMM(kcat_PTP, Km_v24 *omega_cyto); 

fMA(Kf_v25/omega_extra); 

fMA(Kr_v25); 

fMM(kcat_PPase_MAPK, Km_v26 * omega_cyto); 
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kineticLawOf v27 : fMM(kcat_pp2a_4_pp2a_4, Km_pp2a_4 *omega_cyto); 

kineticLawOf v28 : fMM(kcat_pp_ptp_pp_ptp, Km_v28 * omega.cyto); 

kineticLawOf v29 : fMM(kcat_GRK_bg_GRK_bg * cyto_mem /omega_cyto, Km_GRK_bg * cyto_mem); 

kineticLawOf v30f : fMA(Kf G binds JBAR/omega cyto); 

kineticLawOf v30b : fMA(Kr_G_binds_BAR); 

kineticLawOf v31 : fMM(kcat_PPase_mek, Km_v31 * omega.cyto); 

kineticLawOf v32f : fMA(Kf_v32 /omega.cyto); 

kineticLawOf v32b : fMA(Kr_v32); 



AC .active 

G_GDP 

G_protein 

G_a_s 

GRK_bg 

iso_BAR_p 

PDE4 

ATP 

R2C2 

PPJPDE 

BAR 

BAR_G 

iso 

iso_BAR 



MAPK_active = 
MEK 

MEK_active 
B_Raf_active = 
bg 

B_Raf 
PKA 

AC 

AMP 
GRK 



vl6f >> AC .active @cy to _mem + vl6b << AC_active@cyto_mem + 

v09 (+) AC .active @cy to _mem; 

vlOf » G_GDP@cyto + vlOb « G_GDP@cyto + 

vllf « G_GDP@cyto + vllb » G_GDP@cyto; 

vllf » G_protein@cyto + vllb « G_protein@cyto + 

vl2f << G_protein@cyto + vl2b >> G_protein@cyto + 

v30f « G_protein@cyto + v30b » G_protein@cyto; 

vOlf » G_a_s@cyto + vOlb « G_a_s@cyto + vlOf « G_a_s@cyto + 

vlOb » G_a_s@cyto + vl6f « G_a_s@cyto + vl6b » G_a_s@cyto; 

v06f » GRK_bg@cyto + v06b « GRK_bg@cyto + 

v29 (+) GRK_bg@cyto; 

vl9 >> iso_BAR_p@cyto_mem + v29 >> iso_BAR_p@cyto_mem; 
v27 » PDE4@cyto + v20 « PDE4@cyto + v05 (+) PDE4@cyto; 
v09 « ATP@cyto + vl7 « ATP@cyto; 
vl8f « R2C2@cyto + vl8b » R2C2@cyto; 
v27 (+) PPJ>DE@cyto; 

v04f « BAR@cyto_mem + v04b >> BAR@cyto_mem + 
v30f « BAR@cyto_mem + v30b >> BAR@cyto_mem; 
v30f » BAR_G@cyto_mem + v30b « BAR_G@cyto_mem + 
v25f « BAR_G@cyto_mem + v25b » BAR_G@cyto_mem; 
v04f « iso@extra + v04b » iso@extra + v25f « iso@extra + 
v25b » iso@extra; 

vOlf » iso_BAR@cyto_mem + vOlb « iso_BAR@cyto_mem + 
v04f » iso_BAR@cyto_mem + v04b « iso_BAR@cyto_mem + 
vl2f « iso_BAR@cyto_mem + vl2b » iso_BAR@cyto_mem + 
vl9 « iso_BAR@cyto_mem + v29 « iso_BAR@cyto_mem; 
v07 » MAPK_active@cyto + v22 « MAPK_active@cyto + 
v24 « MAPK_active@cyto + v26 « MAPK_active@cyto; 
v31 » MEK@cyto + v21 « MEK@cyto; 
v21 » MEK_active@cyto + v31 « MEK_active@cyto + 
v07 (+) MEK_active@cyto; 

v08 » BJRaf_active@cyto + v03 « B _Raf _active @ cy to + 
v21 (+) B _Raf .active @cy to; 

vOlf » bg@cyto + vOlb « bg@cyto + v06f << bg@cyto + 
v06b » bg@cyto + vl If « bg@cyto + vl lb >> bg@cyto; 
v03 » B_Raf@cyto + v08 « B_Raf@cyto; 

vl3f » PKA@cyto + vl3b « PKA@cyto + v08 (+) PKA@cyto + 

vl5 (+) PKA@cyto + v20 (+) PKA@cyto; 

vl6f « AC@cyto_mem + vl6b >> AC@cyto_mem + 

vl7 (+) AC@cyto_mem; 

v02 » AMP@cyto + v05 » AMP@cyto + vl4 » AMP@cyto; 
v06f « GRK@cyto + v06b » GRK@cyto + vl9 (+) GRK@cyto; 
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PP2A 
MAPK 

PTP 

PTP PKA 

c_R2C2 
c2_R2C2 
c3_R2C2 
iso_BAR_G 



PDE_high_km 
cAMP 



PTP_PP 
PDE4_P 



v03 (+) PP2A@cyto + v26 (+) PP2A@cyto + v31 (+) PP2A@cyto; 
v22 » MAPK@cyto + v24 » MAPK@cyto + v26 » MAPK@cyto + 
v07 « MAPK@cyto; 

v28 » PTP@cyto + vl5 « PTP@cyto + v24 (+) PTP@cyto; 
vl5 » PTP_PKA@cyto + v28 « PTP_PKA@cyto + 
v22 (+) PTP_PKA@cyto; 

vl8f » c_R2C2@cyto + vl8b « c_R2C2@cyto + 

v23f « c_R2C2@cyto + v23b » c_R2C2@cyto; 

v23f » c2_R2C2@cyto + v23b « c2_R2C2@cyto + 

v32f « c2_R2C2@cyto + v32b » c2_R2C2@cyto; 

v32f » c3_R2C2@cyto + v32b « c3_R2C2@cyto + 

vl3f « c3_R2C2@cyto + vl3b » c3_R2C2@cyto; 

vl2f » iso_BAR_G@cyto_mem + vl2b « iso_BAR_G@cyto_mem + 

v25f » iso_BAR_G@cyto_mem + v25b « iso_BAR_G@cyto_mem + 

vOlf « iso_BAR_G@cyto_mem +v01b » iso_BAR_G@cyto_mem; 

vl4 (+) PDE_high_km@cyto; 

v09 » cAMP@cyto + vl7 » cAMP@cyto + v02 « cAMP@cyto + 
v05 « cAMP@cyto + vl3f « cAMP@cyto + vl3b » cAMP@cyto + 
vl4 « cAMP@cyto + vl8f « cAMP@cyto + vl8b » cAMP@cyto + 
v23f « cAMP@cyto + v23b » cAMP@cyto + v32f « cAMP@cyto + 
v32b » cAMP@cyto; 
v28 (+) PTP_PP@cyto; 

v20 » PDE4_P@cyto + v27 « PDE4_P@cyto + v02 (+) PDE4_P; 



MAPK_active Jraction = MAPK_active@cyto / (MAPK@cyto + MAPK_active@cyto); 

MAPK_Pathway :: = B_Raf@cyto[120] < * > B_Raf_active@cyto[0] < * > MEK@cyto[108] < * > 
MEK_active@cyto[0] < * > MAPK@cyto[217] < * > 
MAPK_active@cyto[0] < * > PP2A@cyto[60] < * > PTP@cyto[120] < * > 
PTPJ > KA@cyto[0] < * > PTP_PP@cyto[60]; 

cAMP .Pathway :: = cAMP@cyto[0] < * > ATP@cyto[3010000] < * > PKA@cyto[0] < * > 

PDE4_P@cyto[0] < * > PDE4@cyto[241] < * > PP_PDE@cyto[120] < * > 
PDE_high_km@cyto[301] < * > AMP@cyto[3010000] < * > 
R2C2@cyto[120] < * > c JR2C2@ C yto[0] < * > c2_R2C2@cyto[0] < * > 
c3_R2C2@cyto[0]; 

G_Pathway :: = iso@extra[668] < * > BAR@cyto_mem[19] < * > G_protein@cyto[2167] < * > 

iso_BAR@cyto_mem[0] < * > iso_BAR_p@cyto_mem[0] < * > 
isoJBAR_G@cyto_mem[0] < * > G_a_s@cyto[0] < * > AC@cyto_mem[60] < * > 
AC_active@cyto_mem[0] < * > G_GDP@cyto[0] < * > GRK_bg@cyto[0] < * > 
BAR_G@cyto_mem[0] < * > bg@cyto[0] < * > GRK@cyto[l]; 

G_Pathway <v09, vl7> cAMP .Pathway <v08, vl5> MAPK_Pathway 



